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ABSTRACT 


Solar  Cyclical  Trend  Study  of  the  Mid-latitude, 
Quiet-Time,  Meridional,  Neutral  Winds 
at  Winter  Solstice  Conditions 


Ronald  Lee  Breninger,  Master  of  Science 
Utah  State  University,  1989 


Major  Professor:  Dr.  Kent  L.  Miller 

Department:  Center  for  Atmospheric  and  Space  Sciences 


Located  within  the  region  of  the  thermosphere  is  the 
major  portion  of  the  ionosphere.  Distribution  of  the 
ionospheric  plasma  within  this  region  is  a  function  of 
atmospheric  mass  and  energy  transport.^  At  mid-latitudes, 
absorption  of  UV  and  EUV  radiation  arJ  two  of" the  principle 
ionisation  and  heating  sources  within  the  thermosphere.  The 
level  of  UV  and  EUV  radiation/fs  directly  related  to  the 
amount  and  intensity  of  aeoar  activity.  As  with  temperature 
variations,  other  thermospheric  properties,  such  as  particle 
density,  and  circulation  phenomena,  such  as  zonal  and 
meridional  wind^,  indicate  a  dependence  on  the  level  of 
solar  activity.  Systematic  measurements  of  thermospheric 
winds  are  ma/de  on  a  limited  basis  by  a  sparse  network  of 


Fabry-Peroti  interferometers,  incoherent  scatter  radars,  and 
satellites.^  For  the  current  study,  the  component  of  the  —/ 1 


/-•  f  f  *  *  ^ 
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neutral  wind  along  the  magnetic  meridian  is  derived  from 
ground-based  ionosonde  measurements  of  the  F2  peak  layer 
height.  Meridional  wind  variations  with  respect  to 
location,  universal  time,  and  level  of  solar  activity  are 
the  focal  points  on  this  investigation.  The  primary 
timescale  of  interest  covers  a  period  of  one  solar  cycle, 
from  1977  to  1987.  Data  from  one  station  have  been  extended 
to  1965  to  study  the  repetitive  nature  of  solar  activity  on 
the  meridional  winds. 

The  F2  layer  height  is  determined  by  using  scaled  iono¬ 
sonde  critical  frequencies  and  an  empirical  equation  devel¬ 
oped  by  Bradley  and  Dudeney.  Ionospheric  parameters  relat¬ 
ing  to  the  F2  layer  height  were  determined  from  quiet  to 
moderate  solar  indices  using  the  Field  Line  Interhemispheric 
Plasma  (FLIP)  ionospheric  model  of  Richards  and  Torr. 

^  Results  of  this  study  indicate  a  definite  variation  of 
wind  speed  and  direction,  which  correlates  with  changing 
levels  of  solar  activity.  It  is  hoped  that  the  results  of 
this  study  will  aid  current  efforts  to  develop  ionospheric 
models  and  enhance  their  forecasting  capabilities.  ,  \ 

(81  pages^ 


CHAPTER  I 


INTRODUCTION 


Overview 

For  the  purposes  of  this  study,  the  thermosphere  is  the 
region  above  the  mesopause,  where  the  neutral  gas 
temperature  begins  to  increase  with  height.  It  extends  up 
to  the  exosphere,  where  the  neutral  temperature  is  constant. 
Located  within  the  thermosphere  is  the  major  portion  of  the 
ionosphere,  which,  under  the  influence  of  electric  and 
magnetic  fields,  directly  effects  the  energetics  and 
dynamics  of  the  neutral  atmosphere.  Distribution  of  the 
ionospheric  plasma  within  this  region  is  a  function  of 
atmospheric  mass  and  energy  transport  [Blum  and  Harris, 

1973]  . 

Absorption  of  ultraviolet  (UV)  and  extreme  ultraviolet 
(EUV)  solar  radiation  are  two  of  the  principle  ionization 
and  heating  sources  of  the  mid-latitude  thermosphere.  The 
level  of  UV  and  EUV  radiation  that  reaches  the  earth's  upper 
atmosphere  is  largely  dependent  on  the  amount  and  intensity 
of  solar  activity  [Lakshimi  et  al.,  1988]  and  causes  large 
temperature  variations  in  the  thermosphere,  as  seen  in 
Figure  1  [Banks  and  Kockarts,  1973;  Jacchia  and  Slowey, 

1973;  Hedin  and  Mayr,  1987].  A  similar  solar  cycle  effect 
is  expected  in  other  thermospheric  properties,  such  as 
particle  density,  and  circulation  phenomena,  such  as 


100  SOQ  1000  1900  2000 
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Fig.  l.  Vertical  temperature  profile  of  the  earth's 
atmosphere.  [Reprinted  by  permission  of  Academic 
Press  from  Banks  and  Kockarts  1973]. 


zonal  and  meridional  winds  [Hernandez  and  Roble,  1976; 
Hernandez,  1982]. 

Thermal  energy,  which  is  deposited  in  the  upper  atmos¬ 
phere,  is  transported  by  vertical  heat  conduction  to  the 
lower  atmosphere,  where  radiative  cooling  processes 
dominate.  Forcing  functions  due  to  UV  and  EUV  radiation, 
coupling  pro-cesses  of  the  thermosphere-magnetosphere 
systems,  and  coupling  processes  of  the  thermosphere-lower 
atmosphere  systems  exhibit  spatial  and  temporal  variations. 
A  combination  of  thermospheric  winds  from  solar  radiation 
inputs,  electric  field  convection  processes,  and  ion  drag 
effects  cause  frictional  interactions  of  Joule  heating  and 
viscous  energy  dissipation  [Mayr  and  Harris,  1978;  Mayr  and 
Harris  1979;  Forbes,  1982;  Hedin  and  Mayr,  1987], 

Problem  Statement 

Systematic  measurement  of  thermospheric  winds  is  made 
on  a  limited  basis  by  a  sparse  network  of  Fabry-Perot 
interferometers  [Hernandez  and  Roble,  1976],  incoherent 
scatter  radars  [Salah  and  Holt,  1974],  and  by  satellite 
measurements  of  the  Doppler  shift  of  the  airglow  spectra 
[Herraro  et  al.,  1988],  A  method  of  deriving  the  component 
of  the  neutral  wind  along  the  magnetic  meridian,  using 
ground-based  ionosonde  measurements  of  the  F2  peak  layer 
height  as  developed  by  Miller  et  al.  [1986],  greatly 
enhances  the  capabilities  of  studying  mid-latitude  cir¬ 
culation  phenomena.  Meridional  winds  determined  by  this 
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method  are  a  function  of  location,  universal  time,  and  the 
level  of  solar  activity.  These  variations  are  the  primary 
areas  of  concern  in  this  research  effort. 

Objectives  of  the  Investigation 

The  primary  objective  of  this  investigation  is  to 
conduct  a  trend  study  of  the  mid-latitude,  quiet-time, 
meridional  neutral  winds  over  the  period  of  two  solar 
cycles,  which  span  the  time  period  of  1965  to  1987,  at 
winter  solstice  conditions.  The  use  of  archived  ionosonde 
data,  from  the  World  Data  Center  at  Boulder,  Colorado 
[Conkright  and  Brophy,  1982],  allows  for  such  a  study  to  be 
conducted  on  a  global  scale.  The  method  of  determining 
meridional  neutral  winds  is  based  on  measurement  of  the 
height  of  the  maximum  electron  density  of  the  F  layer 
[Miller  et  al.,  1986].  The  height  of  the  maximum  electron 
density,  hmax  or  hmF2  as  it  is  referred  to,  increases  or 
decreases  as  collisions  between  ions  and  neutral  atmospheric 
particles  force  the  ions  to  move  parallel  to  the  magnetic 
field  in  the  direction  of  the  meridional  component  of  the 
neutral  wind  [Miller  et  al.,  1987].  Natural  downward 
diffusion  of  ionization  is  either  enhanced  or  opposed, 
depending  on  the  direction  of  the  meridional  neutral  wind. 

In  response  to  such  activity,  the  maximum  ion  density  occurs 
at  a  higher  or  lower  level  of  the  ionosphere  with  respect  to 


the  "balance  height”  of  no  neutral  wind. 

A  second  objective  of  this  study  is  to  conduct  a  sta- 
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tistical  error  analysis  of  the  technique  of  Miller  et  al. 
[1986]  for  determination  of  the  meridional  neutral  wind. 

The  error  in  the  meridional  neutral  wind  is  determined  to  be 
a  function  of  time,  change  in  the  layer  height  with  respect 
to  wind  speed,  hmax,  and  the  balance  height,  which  are  the 
variables  used  in  the  calculation  of  the  meridional  winds. 
The  definitions  of  these  variables  are  iterated  in  a  later 
portion  of  this  paper.  A  finite  differencing  technique, 
using  the  root  mean  square  error  of  each  of  the  variables, 
is  used  to  establish  the  uncertainty  in  the  meridional 
neutral  wind,  U  [Baird,  1962]. 

The  third  objective  is  to  develop  a  representation  of 
the  meridional  neutral  winds  by  the  use  of  a  grid  layout  of 
ionosonde  stations.  Stations  selected  for  this  study  form  a 
latitudinal  chain  and  a  longitudinal  chain,  as  seen  in 
Figure  2.  Additional  site-specific  information  is  contained 
in  Appendix  A.  Such  a  network  of  stations  allows  for  the 
investigation  of  the  meridional  neutral  winds  in  a  spatial 
coordinate  frame  [Challinor  and  Eccles,  1971].  The  result 
of  this  study  provides  location-specific  variations  of  the 
meridional  neutral  winds  with  respect  to  solar  cycle 
variations . 


6 


e 

o 

u 

<w 


<D  TJ 
P  <V 
«J  P 

c  a 
•P  o 

TJ  TJ 
p  «J 

o  •— < 
o 

o  c 
o 

Q)  -H 
i — I  p 

o  o 

D,  0) 

•H  -l—l 


TJ  O 

p 

•a  a 


<u 

p 

a> 


p  p 

C  <0 

d  o 
o  a) 
e 

(0 

o 

C  -p 

o  x: 


a 
(0 
p 

O' 
•H  o 
p  a) 

<o  tP- 


oi 

c 

o 


o 

o 


in 

(0  VD 
CTv 

C  «-l 

o 


0> 

TJ 
C  TJ  • 
O  <DH 

tn  w  «j 
o  o 

Q.4J 
P  Q> 
Q> 

a  >, 

3  M 
0)  (0 
o 

fl>  C 

•  e  p 

O'  <0  ® 
•P  PJ3 
Cn  P  O 


c 

o 


OJ 


7 


CHAPTER  II 

REVIEW  OF  LITERATURE 


Meridional  Winds 

Thermal  expansion  of  the  atmosphere  by  day  forms  a 
diurnal  bulge  [Rishbeth,  1972;  Salah  and  Holt,  1974],  which 
gives  rise  to  horizontal  pressure  gradients,  which  in  turn 
drive  horizontal  winds.  These  thermospheric  winds,  to  some 
degree,  are  subject  to  the  Coriolis  effect  due  to  the 
earth's  rotation.  They  are  even  more  strongly  influenced  by 
frictional  forces  due  to  the  viscosity  of  the  air  and  col¬ 
lisions  between  ions  and  neutral  atmospheric  particles. 
Additional  information  on  this  subject  is  provided  in  a 
later  portion  of  this  chapter.  Ions  exert  a  drag  on  the  air 
motion  because  of  the  strong  ties  with  the  earth's  magnetic 
field.  Ion  drag  is  a  major  limiting  factor  of  the  wind 
speeds  at  thermospheric  heights  [Rishbeth,  1972]. 

Thermospheric  winds  blow  from  the  hottest  part  of  the 
thermosphere  to  the  coldest  part  and  therefore  blow  across 
the  poles  and  zonally  around  the  earth  in  the  low  latitude 
regions.  The  winds  move  F- region  ions  and  electrons  in  the 
direction  of  the  earth's  magnetic  field,  with  the  overall 
effect  of  the  wind  being  dependent  on  its  orientation  with 
respect  to  the  magnetic  field. 

Daytime  (poleward)  and  nighttime  (equatorward)  meridi¬ 
onal  winds  cause  a  lowering  of  the  daytime  F2  layer  and 
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raising  of  the  nighttime  F2  layer  of  the  ionosphere  beyond 
the  normal,  diurnal  variations  caused  by  production,  recom¬ 
bination,  and  diffusion  processes  [Rishbeth,  1966;  Rishbeth, 
1972].  The  variation  of  the  F2  peak  layer  height  due  to  the 
influence  of  the  meridional  neutral  winds  is  nearly  linear, 
under  steady  state  conditions  [Miller  et  al.,  1986;  Miller 
and  Torr,  1987].  The  neutral  wind,  as  determined  by  such  a 
method,  is  given  by  the  relationship 


where  U  *  horizontal  component  of  the  neutral  wind  along  a 
magnetic  meridian, 

hmax  s  measured  F2  peak  layer  height, 

ho  s  calculated  balance  height  of  the  F2  peak  density 
at  U  =  0,  as  calculated  by  the  Richards-Torr 
ionospheric  model  [Young  et  al.,  1980], 


a  3  calculated  value  of  (hHI2  -  hKT1)/(u2  -  Uj) ,  as 

determined  by  a  linear  regression  technique  of  the 
Richards-Torr  ionospheric  model  [Young  et  al., 
1980] , 

hRT1,  hKT2  b  layer  height  as  determined  by  the  Richards- 
Torr  ionospheric  model  [Young  et  al., 

1980],  and 

ulf  u2  ■  inputed  meridional  winds  of  the  Richards- 

Torr  ionospheric  model  [Young  et  al.,  1980]. 

Derivation  of  equation  (l)  neglects  electric  field  effects 
on  the  F2  layer  height.  There  is  no  component  of  the  elect¬ 
ric  field  parallel  to  the  magnetic  field  at  mid-latitudes 
which  affects  the  layer  height  [Stubbe  and  Chandra,  1970; 
Rishbeth  1972;  Miller  et  al.,  1986;  Miller  et  al.,  1987; 
Miller  and  Torr,  1987].  At  mid-latitudes,  it  is  the  east- 
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west  component  of  the  earth's  electric  field  which  affects 
the  F2  layer  height.  The  apparent  meridional  wind,  U,  is 
therefore  a  combination  of  the  meridional  wind,  Uy,  and  the 
zonal  electric  field  wind  [Hanson  and  Patterson,  1964; 
Hernandez  and  Roble,  1976;  Hiller  et  al.,  1986;  Miller  et 
al.,  1987],  or 

u  =  Uy  +  E,  (2) 

B  sin  I 

For  the  northern  hemisphere,  Uy  is  the  magnetically  north¬ 
ward  component  of  the  actual  neutral  wind;  Ez  the  eastward 
component  of  the  electric  field,  E;  B  the  magnitude  of  the 
magnetic  field,  B;  and  I  the  inclination  of  the  magnetic 
field.  At  mid-latitudes,  the  effect  of  the  second  term  of 
equation  (2)  is  generally  small  in  comparison  to  the  first 
term  [Stubbe  and  Chandra,  1970;  Salah  and  Holt,  1974;  Miller 
et  al.,  1987].  Additional  details  on  this  point  are  includ¬ 
ed  in  the  section  on  ionospheric  models.  The  area  of  great¬ 
est  uncertainty  of  U  occurs  at  sunrise  due  to  the  effects  of 
photochemistry  on  the  F2  layer  height  [Hanson  and  Patterson, 
1964;  Rishbeth,  1966].  Photoionization  establishes  a  mini¬ 
mum  value  of  hmax  at  this  time  of  the  day  [Miller  et  al., 
1986].  The  daytime  F2  layer  height  is  most  sensitive  to 
changes  in  the  meridional  wind,  since  there  is  a  balance 
between  diffusive  and  photochemical  processes  occurring 
within  the  ionosphere.  It  is  less  sensitive  at  night,  when 
the  F  layer  is  maintained  by  downward  diffusion. 


The  major  contributor  to  thermospheric  circulation,  at 
mid-latitudes,  is  heating  by  absorption  of  solar  UV  and  EUV 
radiation  [Rishbeth,  1972;  Dickinson  et  al.,  1977;  Hedin  and 
Mayr,  1987],  with  some  amplification  resulting  from  high 
latitude  heat  sources.  At  winter  solstice  conditions,  a 
summer  to  winter  interhemispheric  flow  regime  is  estab¬ 
lished.  This  flow  regime  is  reinforced  by  a  high  latitude 
heat  source  in  the  summer  hemisphere  and  opposed  in  the 
winter  hemisphere  at  F  region  heights  above  150  km  [Roble  et 
al.,  1977;  Babcock  and  Evans,  1979].  Below  150  km  a  summer 
to  winter  flow  pattern  is  established  at  all  latitudes. 

Solar  UV  and  EUV  heating  is  approximately  two  times 
lower  at  solar  minimum  than  it  is  at  solar  maximum.  Cor¬ 
responding  to  this  solar  cycle  effect,  a  factor  of  4.5 
difference  has  been  noted  in  the  strength  of  the  high  lati¬ 
tude  heat  source  between  solar  maximum  and  solar  minimum 
conditions  [Roble  et  al.,  1977].  This  variation  is  computed 
based  on  a  requirement  for  modeling  of  the  temperatures  and 
meridional  winds  to  match  observations  made  at  solar  cycle 
maximum  and  minimum  conditions.  The  high  latitude  heat 
source  variation  is  the  result  of  the  interaction  on  the 
solar  wind  and  the  earth's  magnetic  field  [Dickinson  et  al., 
1977] . 

In  a  study  by  Roble  et  al.  [1977],  at  winter  solstice, 
a  temperature  difference  of  325  K,  at  solar  maximum,  and 
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200  K,  at  solar  minimum,  occurs  between  the  summer  and 
winter  poles.  Latitudinal  variations  of  the  mean  molecular 
mass  are  also  greatest  at  solstice  conditions  with  a  maximum 
value  occurring  in  the  summer  hemisphere.  This  variation  is 
due,  in  part,  to  the  large  mixing  ratios  of  the  heavier  at¬ 
mospheric  constituents  found  in  the  winter  hemisphere. 

Under  solstice  conditions,  at  solar  maximum,  the  entire 
thermospheric  circulation  below  150  km  exhibits  the  summer 
to  winter  hemispheric  flow  pattern  as  seen  in  Figure  3. 

Above  200  km,  the  meridional  wind  strength  intensifies  in  an 
eguatorial  direction  with  a  maximum  strength  occurring  at 
the  solstice  time  period  [Roble  et  al.,  1977].  Under  sol¬ 
stice  conditions,  the  auroral  zone  heating  in  the  summer 
hemisphere  combines  with  solar  EUV  heating  to  drive  a  Hadley 
cell  circulation,  which  transports  air  across  the  equator 
from  the  summer  to  winter  hemisphere.  This  flow  is  resisted 
by  an  oppositely  directed  flow  in  the  winter  hemisphere  due 
to  auroral  heating  in  that  hemisphere  [Babcock  and  Evans, 
1979].  The  reverse  cell  acts  to  decrease  the  poleward  winds 
in  the  winter  hemisphere,  which  are  driven  by  solar  heating 
in  the  summer  hemisphere,  and  results  in  an  assymetric 
pattern  in  the  meridional  winds  at  mid-latitudes  [Roble  et 
al.,  1977;  Babcock  and  Evans,  1979].  The  strength  of  this 
circulation  and  the  location  of  the  boundary  of  circulation 
reversal,  as  illustrated  in  Figure  3,  are  expected  to  be 
dependent  on  the  level  of  geomagnetic  activity  [Hernandez 
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and  Roble,  1977;  Roble  et  al.,  1977], 

Thermospheric,  meridional  winds  at  solar  minimum  and  at 
winter  solstice  conditions  flow  from  the  summer  to  winter 
hemisphere  at  all  altitudes  and  latitudes  [Hernandez  and 
Roble,  1977].  The  magnitude  of  these  winds,  as  seen  in 
Figure  4,  is  nearly  two  times  lower  than  the  values  found  at 
solar  maximum.  This  difference  is  attributable  to  the 
diminished  high  latitude  heat  source,  which  is  needed  to 
create  a  reverse  circulation  [Roble  et  al.,  1977]. 

Observations  of  sunspot  numbers  provide  a  means  of 
prediction  of  ionospheric  parameters,  such  as  the  F  region 
critical  frequency,  foF2,  and  the  F  region  peak  layer 
height,  hmF2  [Smith  and  King,  1981].  Drawbacks  to  using  the 
monthly  mean  sunspot  number,  R12/  are  the  saturation  effects 
and  the  variation  of  foF2  for  the  same  value  of  R12  in 
different  phases  of  the  solar  cycle  [Hedin  and  Mayr,  1987; 
Lakshimi  et  al.,  1988].  Results  of  a  study  by  Lakshimi  et 
al.  [1988]  of  the  values  of  foF2  versus  EUV  flux  (170-190  A) 
and  the  sunspot  number,  RIZ,  show  no  saturation  effects  in 
the  winter  solstice  time  frame  at  low  latitudes.  Such  is 
not  the  case  for  mid-latitude  regions,  as  neither  the  use  of 
EUV  radiation  flux  or  R12  provide  a  greater  degree  of 
accuracy  in  the  determination  of  foF2.  Additionally,  the 
results  of  the  work  of  Babcock  and  Evans  [1979]  and  Lakshimi 
et  al.  [1988]  indicate  that  an  increase  in  solar  activity 
increases  the  pressure  gradients,  which  drive  thermospheric 


Fig.  4.  Contours  of  temperature  and  winds  for  winter 
solstice  at  solar  minimum;  (a)  meridional  winds  (m/s) 
(b)  zonal  winds  (m/s) ,  (c)  temperature  (K) ,  and  (d) 
vertical  winds  (cm/s)  [from  Roble  et  al . ,  1977]. 
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winds,  but,  at  the  same  time,  an  increase  in  the  charged 
particle  population  leads  to  increased  ion  drag,  which 
inhibits  the  winds  [Hedin  and  Mayr,  1987].  In  general,  as 
the  level  of  solar  activity  decreases,  daytime  wind  speeds 
increase  because  ion  drag  decreases  more  rapidly  than  the 
pressure  gradient  force  [Rishbeth,  1972], 

Based  on  a  total  of  64  geomagnetical ly  quiet  days,  from 
1970  to  1975,  Babcock  and  Evans  [1979]  found  a  decrease  in 
the  annual  mean  meridional  wind,  at  Millstone  Hill,  of  25 
m/s  equatorward  to  0  m/s.  They  conclude  that  the  decrease 
in  wind  speed  is  indicative  of  a  decrease  in  auroral  forcing 
relative  to  the  level  of  solar  forcing  at  solar  minimum. 
Babcock  and  Evans  [1979]  carried  out  a  least  squares  curve 
fitting  of  the  mean  meridional  wind,  at  the  Millstone  Hill 
observatory  in  Boston,  Massachusetts,  and  established  a 
correlation  with  the  three  month  average  10.7  cm  solar  radio 
flux,  F10.7A  [Jacchia  and  Slowey,  1973;  Hedin  and  Mayr, 
1987].  At  Millstone  Hill,  for  solar  maximum  the  mean  summer 
wind  is  approximately  75  m/s  equatorward  while  a  winter  wind 
of  30  m/s  poleward  was  calculated.  For  solar  minimum,  the 
average  winter  wind  increases  to  40  m/s  and  the  duration  of 
the  summer  equatorward  wind  is  reduced  in  comparison  to  that 
of  solar  maximum  conditions  while  the  speed  was  relatively 
unchanged . 

Ionospheric  Models 

In  the  thermosphere,  molecules  collide  frequently  so 
that  the  air  may  be  considered  to  be  a  fluid,  which  is  sub- 
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ject  to  the  hydrodynamic  equations  of  motion.  Furthermore, 
the  air  can  be  treated  as  a  single  fluid  because  the  dif¬ 
ferential  motion  of  the  constituents  is  much  less  than  the 
overall  wind  speed  [Rishbeth,  1972;  Blum  and  Harris,  1973]. 

The  equation  of  motion  may  be  represented  as: 

dD  +  2n  x  u  =  7P  -  i/ni(U  -  V)  +  it  v2U  +  g  (3) 

dt  p 

(a)  (b)  (c)  (d)  (e)  (f) 

where  term  (a)  is  acceleration,  (b)  Coriolis  term,  (c) 

pressure  gradient,  (d)  ion  drag,  (e)  viscous  drag,  and  (f) 

gravity.  For  the  terms  of  equation  (3)  U  is  the  neutral 

wind  velocity,  V  is  the  ion  drift  velocity,  n  the  earth's 

angular  velocity,  uni  the  neutral-ion  collision  frequency, 

H/p  the  kinematic  viscosity,  and  g  the  acceleration  due  to 

gravity. 

In  order  to  calculate  U(t)  in  a  coordinate  frame  which 
is  fixed  with  respect  to  the  earth,  the  Coriolis  and  centri¬ 
petal  accelerations  due  to  rotation  of  the  earth  must  be 
included  in  the  equation  of  motion  for  determination  of  the 
motion  of  the  air,  dU/dt.  Local  acceleration,  au/at, 
differs  from  this  quantity  due  to  the  transport  of  momentum 
past  a  given  point.  According  to  the  chain  rule  for 
derivatives,  d/dt  =  a/at  +  (U  •  V) .  A  non-linear  advection 
term,  (U  •  v)U,  results  from  application  of  this  definition 
to  the  wind  velocity  U. 

Centripetal  acceleration  is  generally  neglected  because 
it  is  constant  in  magnitude  and  direction  at  any  given 
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latitude.  The  non-linear  advection  term  is  small  if  U  is 
much  smaller  than  the  earth's  rotational  speed,  n,  and 
therefore  may  generally  be  neglected.  The  non-linear 
advection  term  may  be  significant  for  large  values  of  U  and 
large  spatial  gradients,  V  •  U,  as  can  occur  near  the 
periods  of  sunrise  and  sunset  [Rishbeth,  1972;  Salah  and 
Holt,  1974]. 

The  pressure  gradient  force  is  not  exactly  vertical  and 
therefore  is  comprised  of  both  horizontal  and  vertical  com¬ 
ponents.  It  is  the  horizontal  components  which  provide  the 
driving  force  for  zonal  and  meridional  winds.  The  vertical 
component,  though  largest  in  magnitude,  is  counterbalanced 
by  the  gravitational  force  and  has  little  effect  on  air 
motion . 

Viscous  forces  are  created  if  wind  shears  exist.  It  is 
generally  the  vertical  wind  shear  which  is  of  concern  and 
therefore  the  viscous  drag  term  of  equation  (3)  becomes 
(M/ p)8zJJ/azz ,  which  allows  for  diffusion  effects  to  be 
considered.  In  other  words,  the  momentum  of  air  tends  to 
diffuse  so  as  to  smooth  the  velocity-height  profile.  This 
term  is  important  at  high  altitudes  where  the  density  of  the 
atmospheric  particles  is  small  and  the  kinematic  viscosity 
is  therefore  large  [Blum  and  Harris,  1973;  Salah  and  Holt, 
1974;  Hedin  and  Mayr,  1987]. 

The  collision  frequency  for  momentum  transfer,  i/ni,  is 
comparable  to  the  average  frequency  with  which  a  neutral 
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particle  experiences  a  collision  with  an  ion.  By  Newton's 
Third  Law,  pni/nl  =  where  />*  and  pn  are  ion  and  neutral 

densities.  At  the  F2  peak  altitude,  where  1  part  in  1000  of 
the  air  is  ionized,  uiB  »  1  s'1  and  uBi  «  10'3  s'1.  The  time 
constant  within  which  an  ion  transfers  motion  to  the  neutral 
is  l/vnl  »  1000  s  [Rishbeth,  1972]. 

The  four  principal  causes  of  ion  motion  in  the  F  layer 
are  plasma  diffusion,  neutral  horizontal  winds,  vertical 
thermal  expansion  and  contraction  of  the  atmosphere,  and 
drifts  due  to  electric  fields  generated  by  dynamo  action  in 
the  E  layer.  Their  relative  magnitudes,  in  terms  of  their 
contributions  to  V,  at  mid-latitudes  are  given  in  Table  1 
[Rishbeth  1972;  Salah  and  Holt,  1974].  If  the  ion  gyrofre- 
quency,  wlt  is  much  greater  than  the  collision  frequency, 
i/in,  as  is  the  case  at  F  layer  altitudes,  then  plasma  dif¬ 
fusion,  the  neutral  horizontal  winds,  and  thermal  expansion 
and  contraction  produce  ion  motion  which  is  parallel  to  the 
geomagnetic  field  [Roble  and  Dickinson,  1974].  Ion  drift 
due  to  the  electric  field  is  perpendicular  to  the  magnetic 
field,  B. 

An  approximation  of  the  ion  motion  is  made  by  setting 
the  ion  velocity  equal  to  the  meridional  component  of  the 
neutral  wind,  U,  or 

V  =  (P  »  BIB  (4) 

Bz 

Ion  drag  then  depends  only  on  the  orientation  of  the  mag¬ 
netic  field  and  is  a  function  of  the  inclination,  I,  and 


Table  1.  Magnitude  of  the  Principle  Causes  of  Ion 
Motion  Relative  to  the  Ion  Drift  Velocity  [based 
on  Rishbeth,  1972]. 
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CAUSE  OF  ION  MOTION 

Plasma  Diffusion 

Neutral  Horizontal 
Minds 

Vertical  Thermal 
Expansion  and 
Contraction 

Dynamo  Induced 
Drifts 


MAGNITUDE 
10  m/s 

100  m/s  (nighttime) 
30  m/s  (daytime) 

3  m/s 
30  m/s 


declination,  D,  angles  [Rishbeth,  1972;  Roble  and  Dickinson, 
1974] . 

Plasma  diffusion  must  be  accounted  for  due  to  electro¬ 
static  force  between  ions  and  electrons.  Ion  diffusion 
velocities  are  the  sum  of  the  ion  and  electron  partial 
pressures.  Collisions  between  electrons  and  ions  and 
between  electrons  and  neutrals  have  little  effect  on  the  ion 
drift  velocity  and  the  overall  air  motion.  This  is  due  to 
the  large  differences  in  mass  of  the  electrons  as  compared 
to  that  of  the  ions  and  neutral  particles. 

Based  on  the  above  discussion,  a  modified  equation  of 
motion  can  be  given  as 

dE  *  VP  +  a  ifu  -  i/ni(U  -  V)  -  2fl  x  U  +  g  (5) 

dt  p  8z* 

which  can  be  resolved  into  horizontal  and  vertical  com- 
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ponents.  The  vertical  component  reduces  to  the  hydrostatic 
equation  and  is  generally  neglected  in  ionospheric  models 
[Rishbeth,  1972;  Young  et  al.,  1980].  In  such  cases  a  model 
atmosphere,  such  as  MSIS  [Hedin,  1983],  is  adopted  so  that 
density  is  treated  as  a  known  quantity.  The  Coriolis 
effect,  relative  to  the  ion  drag  effect,  is  not  a  major 
factor  at  night,  when  the  peak  density,  NmF2,  becomes  small. 
In  general,  when  f  =  2SI  sin  <p,  in  which  <p  is  the  geographic 
latitude,  is  greater  than  i/ni  then  the  Coriolis  effect  is  an 
important  factor.  This  situation  is  found  below  150  km  and 
above  400  km.  However,  above  400  km,  viscosity  outweighs 
both  the  Coriolis  effect  and  ion  drag  [Rishbeth,  1972]. 

The  Richards-Torr  ionospheric  model  makes  use  of 
modified  momentum  and  energy  equations  [Hanson  and  Patter¬ 
son,  1964;  Rishbeth,  1972;  Schunk,  1974;  Hernandez  and 
Roble,  1976;  St.  Maurice  and  Schunk,  1977;  Schunk  and  Nagy, 
1978],  Such  modifications  restrict  the  use  of  this  model  to 
the  mid-latitude  ionosphere,  which  is  a  region  of  closed 
magnetic  field  lines.  Additionally,  the  mid-latitudes  are  a 
region  where  electric  field  convection  effects  are  small  in 
comparison  to  that  of  meridional  winds  [Stubbe  and  Chandra, 
1970].  For  this  model,  the  plasma  consists  of  the  two  major 
ions  of  0+  and  H+,  electrons,  and  some  minor  ions  in  a 
neutral  atmosphere  whose  major  constituents  are  0,  02,  N2, 
and  H.  Modifications  of  the  momentum  and  energy  equations 
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are  based  on  the  following  conditions  [Young  et  al.,  1980]: 

1)  Species  temperature  and  flow  velocity  differences 
are  small,  allowing  one  to  neglect  the  stress  and  non¬ 
linear  acceleration  terms  and  to  use  linear  collision  terms. 

2)  Density  and  temperature  gradients,  perpendicular  to 
the  earth's  magnetic  field,  are  neglected. 

3)  Ion  and  electron  temperature  distributions  are 
isotropic. 

4)  Electric  field  convection  of  the  flux  tube  is 
omitted  and  restricts  useage  to  quiet  or  moderate  geomag¬ 
netic  conditions  [Miller  et  al.,  1987].  Such  conditions  are 
maintained  in  the  current  investigation,  as  seen  in 
Appendix  B. 

The  final  result  is  a  set  of  six  equations  consisting 
of  a  photoelectron  two  stream  equation,  ion  and  electron 
energy  equations,  ion  and  electron  momentum  equations,  and 
the  continuity  equation  [Young  et  al.,  1980],  The  Richards- 
Torr  ionospheric  model  simulates  interhemispheric  plasma 
flow  in  an  entire  flux  tube,  which  spans  the  mid-latitude 
between  magnetically  conjugate  points,  in  the  F  region  of 
the  ionosphere,  as  illustrated  in  Figure  5  [Young  et  al., 
1980].  Of  an  additional  note,  neutral  atmospheric 
parameters  within  the  model  are  calculated  by  the  MSIS-83 
model  [Hedin,  1983].  Selection  of  the  Richards-Torr 
ionospheric  model  is  based  on  its  incorporation  of  seasonal, 
solar  cycle,  and  diurnal  variations.  An  additional 
criterion  is  its  ability  to  respond  to  varying  intensities 


1»  1*:  Boundary  region  (Density  is  continuous) 

2,  2*:  Boundary  region  (Density  and  flux  are  continuous) 

A,  A*:  Regions  of  local  chemistry 

B,  B*:  Regions  of  ion  diffusion  (parabolic  density  equations) 


C  :  Region  of  diffusive  equilibrium 


Fig.  5.  Regions  and  boundaries  of  the  Richards-Torr 
ionospheric  model.  [Reprinted  by  permission  of  Pergamon 
Press  from  Young  et  al . ,  1980]. 
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of  geomagnetic  disturbances  [Young  et  al . ,  1980]. 

The  Bradley-Dudeney  model  ionosphere  provides  a  means 
of  determining  hmax,  or  hmF2,  from  ionosonde  measurements  of 
the  ionosphere.  The  three  parameters,  which  are  inputs  into 
this  model,  are  the  critical  frequencies,  foE  and  foF2,  and 
the  maximum  usable  frequency  factor  for  a  path  of  3000  km 
for  transmission  by  the  F2  layer  of  the  ionosphere, 
M(3000)F2.  Measurement  of  these  parameters  is  made  on  a 
continuous  24-hour  basis  by  over  126  stations,  which  store 
their  data  in  the  World  Data  Center  at  Boulder,  Colorado. 

Figure  6  provides  a  schematic  of  the  Bradley-Dudeney 
model  ionosphere  [Bradley  and  Dudeney,  1973]  which  consists 
of  parabolic  E  and  F  layers  with  a  linear  region  of  electron 
concentration  between  them.  The  calculation  of  hmax  from 
this  model  is  restricted  to  values  of  x,  the  ratio  of  foF2 
to  foE,  greater  than  1.7.  In  this  model,  hmax  is  defined  as 

hmax  =  a  [M(3000)F2]b  (6) 


where 

a  =  1890  -  355 

X  -  1.4 

and 

b  =  (2 . 5x  -  3) "z'35  -  1.6 

Discrepancies  with  the  true  height  profiles  are  of  the  order 
of  ±  10  km  [Bradley  and  Dudeney,  1973;  Dudeney,  1983]  for 
the  quiet,  daytime,  mid-latitude  ionosphere.  Application  of 
this  model  is  based  on  the  selection  of  such  locations  and 
geomagnetic  conditions.  The  error  in  the  Bradley-Dudeney 
formula  of  equation  (6)  is  likely  to  be  greater  than 


j 


Plasm  frequency 


Fig.  6.  Bradley-Dudeney  ionosphere.  The 
ordinate  and  abscissa  are  linear  in  height  and 
frequency.  [Reprinted  by  permission  of  Pergamon 
Press  from  Bradley  and  Dudeney ,  1973]. 
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±  10  km  for  nighttime  conditions  due  to  its  dependency  on 
the  ratio  of  foF2  to  foE.  Nighttime  measurements  of  foE  are 
not  available  due  to  the  decay  of  the  E  region.  For  this 
and  other  problems  of  measurement  of  foE,  such  as  sporadic 
E,  the  values  of  foE  used  in  this  investigation  are 
approximated  from  the  IRI  model  densities  as 

foE  -  9  x  10'3  J  NmE  (7) 

where  NmE  is  the  electron  number  density  at  the  E  peak  layer 
height  [Chen,  1984]. 
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CHAPTER  III 

METHOD 


Procedure 

As  stated  in  the  title,  the  primary  objective  of  this 
investigation  is  a  trend  study  of  the  mid-latitude,  quiet¬ 
time,  meridional  neutral  winds.  The  period  of  study  covers 
two  solar  cycles,  during  the  years  of  1965  to  1987,  with  a 
majority  of  the  emphasis  on  the  years  of  1977  to  1987.  Use 
of  the  additional  solar  cycle,  1965  to  1976,  will  be  as 
verification  of  information  gleaned  from  the  primary  period 
of  this  study.  In  terms  of  seasonal  effects  on  the  thermo¬ 
spheric  winds,  only  the  winter  solstice  time  period,  and  in 
particular  the  month  of  January,  will  be  considered. 

The  method  of  determining  the  values  of  the  meridional 
wind  is  based  on  the  works  of  Rishbeth  [1972]  and  Miller  et 
al.  [1986],  Additional  information  on  this  subject  is 
contained  in  Chapter  II  in  the  section  on  meridional  winds. 
Figure  7  provides  a  schematic  flow  digram  of  the  method 
employed  in  the  determination  of  the  meridional  neutral 
winds . 

As  stated  in  equation  (1),  a  linear  relationship  exists 
between  the  meridional  neutral  wind,  U,  and  the  F2  peak 
layer  height,  hmax  or  hmF2.  Determination  of  the  height, 
hmax,  is  based  on  the  empirical  relationship,  developed  by 
Bradley  and  Dudeney  [1973],  as  stated  in  equation  (6). 


Ionosonde  Data 
(foE,  foF2,  M(3000)F2) 


IRI  - 

(Calculated  foE) 


Bradley  -  Dudeney  Model 
*  hmax 

*  (IRI  hmax  values  are  used  for  years  in  which  no  iono 
sonde  data  is  available.)  I 


Richards  -  Torr  Ionospheric  Model 
a  and  ho 


1 


Miller  Meridional  Hind  Model 
U 


1 

Solar  Activity  Comparison 


Fig.  7.  Schematic  of  data  analysis  for  the  solar  cy¬ 
clical  trend  study  of  the  meridional  neutral  winds. 
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Archived  ionosonde  data,  from  the  World  Data  Center  at 
Boulder,  Colorado,  of  the  critical  frequencies,  foE  and 
foF2 ,  and  the  maximum  transmission  frequency,  M(3000)F2, 
provide  the  observed  input  for  this  investigation.  However, 
due  to  problems  in  the  measurement  of  foE,  as  previously 
stated,  the  E  layer  critical  frequency  used  throughout  this 
investigation  is  calculated  from  the  IRI  model  E  layer 
density,  NmE.  The  use  of  such  data  is  based  on  comparison 
of  the  calculated  values  and  measured  values  of  foE  in  terms 
of  their  relative  impact  on  the  ratio  of  foF2  to  foE,  which 
is  required  for  determination  of  hmax.  The  variance  between 
the  IRI  model  values  and  measured  values  of  foE  is  generally 
less  than  5  to  10%  and  does  not  significantly  change  the 
final  wind  values.  For  the  locations  and  time  periods  where 
ionosonde  data  are  unavailable,  as  stated  in  Appendix  A,  the 
IRI  model  value  of  hmax  has  been  used  in  the  determination 
of  the  meridional  winds  [Miller  et  al.,  1989],  Of  and 
additional  note,  throughout  this  investigation  an  average  of 
the  ten  most  quiet  days  is  used  for  the  determination  of 
hmax,  ho,  and  a.  This  technique  applies  to  the  calculation 
of  both  the  input  ionosonde  frequencies  of  the  Bradley- 
Dudeney  model  and  the  geomagnetic  indices  used  by  the 
Richards-Torr  ionospheric  model.  Selection  of  the  ten  most 
quiet  days  show  no  significant  variance  in  terms  of  Kp 
index,  as  seen  in  Figure  8. 


YEARS 


Fig.  8.  Ten  day  average  of  January  quiet  time  Kp 
index.  Indices  show  no  solar  cycle  variation. 
Error  bars  represent  the  standard  deviation  about 
the  mean. 


30 

As  discussed  in  Chapter  II,  the  Richards-Torr  iono¬ 
spheric  model  is  used  in  this  investigation  for  determina¬ 
tion  of  the  variables  a  and  ho,  which  are  required  for 
calculation  of  the  meridional  winds.  The  development  of  the 
Richards-Torr  ionospheric  model  restricts  its  applicability 
to  the  mid-latitudes,  which  is  the  geographic  region  of 
concern  for  this  investigation.  As  previously  stated,  this 
model  incorporates  seasonal,  solar  cycle,  and  diurnal  varia¬ 
tions.  All  of  these  variations  come  into  play  in  the  de¬ 
termination  of  the  meridional  neutral  wind  [Roble  et  al., 
1977;  Young  et  al.,  1980].  Additionally,  the  Richards-Torr 
model  is  able  to  handle  varying  levels  of  geomagnetic 
activity.  This  is  not  a  critical  criteria  in  the  present 
study  as  only  quiet  to  moderate  levels  of  geomagnetic  ac¬ 
tivity  (Ap  <  8)  are  encountered. 

By  the  use  of  a  grid  system,  which  consists  of  a  long¬ 
itudinal  and  latitudinal  chain  of  ionosonde  stations,  as 
indicated  in  Figure  2,  variations  of  the  meridional  wind 
with  respect  to  location,  universal  time,  and  local  time  are 
studied.  The  timescale  of  this  study  provides  a  means  of 
investigating  the  transition  from  solar  maximum  winds  to 
solar  minimum  winds. 

Data  Analysis 

As  with  most  studies  of  this  magnitude,  two  problem 
areas  must  be  addressed.  First,  a  reduction  of  the  final 
results  to  a  managable,  but  representative,  level  is  es- 
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sential  so  as  not  to  obscure  their  meaning.  Second,  a 
degree  of  statistical  confidence  is  also  required  for  the 
application  of  the  results  to  be  of  use  in  future  studies. 
Only  the  second  problem  area  will  be  addressed  at  the 
present  time  as  the  manner  of  presentation  of  results  will 
be  dealt  with  in  conjunction  with  the  discussion  of  the 
results  of  this  investigation. 

Based  on  equation  (1) ,  one  can  establish  the  error,  or 
uncertainty,  in  the  calculated  meridional  winds  as  a  func¬ 
tion  of  time,  a,  ho,  and  hmax.  By  finite  difference  tech¬ 
niques,  using  the  root  mean  square  method,  the  uncertainty 
in  U  is  given  as: 

-  hmaxV  (6a)2 

“  ) 

or 

SU  -  1 _  [(tfhrnax)2  +  (6ho)2  +  /ho  -  hmaxV  ( 6a) 2  I  1/2 

U  ho  -  hmax  I  ^  a  )  (9) 

Determination  of  6ho  and  6a  are  based  on  the  evaluation  of 

the  method  of  calculation  of  a  and  ho  by  the  Richards-Torr 

ionospheric  model.  For  the  value  of  6 hmax,  from  the 

Bradley-Dudeney  model,  a  value  of  ±  10  km  is  used  in  the 

present  investigation  [Bradley  and  Dudeney,  1973;  Dudeney, 

1983].  By  a  linear  regression  of  a,  as  illustrated  in 


6U 


=  j~  1  |~f  6  hmax) 2  +  (6ho)2  +  ^ho 


L 


+e 


Figure  9, 


By  definition  of  ho. 


■Q  ~  Mi  (^Rt2  “  hm) 


(10) 


ho  *  hRT1  + 


(11) 
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then 

<Sho  = 


/  u2  V  (5hm)2  +  /  ut  \2  (*hRT2)2l  1/2 

( u2  -  uxy  fu2  -  uj 


(12) 


In  accordance  with  the  model  execution  notes  £hm 
$hRT  =  ±  10  km  and  by  substitution 
Sho  =  ( 5hRj) 


—  <5hRT2  — 


"  V  +  Ui 
(u2  -  u,)2 


1/2 


(13) 


The  above  error  analysis  is  strictly  a  statistical  error 
analysis  and  does  not  account  for  the  fact  that  the  model 
is,  at  best,  an  approximation  to  average  conditions. 


1 


^VELOCITY  (m/s) 


Pig.  9.  Statistical  error  analysis  of  a  and 
ho  by  linear  reression  techniques  of  the 
Richards-Torr  ionospheric  model . 
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CHAPTER  IV 

RESULTS  AND  DISCUSSION 

Solar  Cvcle  Variations 

The  major  contributor  to  the  mid-latitude  circulation 
pattern  is  heating  by  absorption  of  solar  UV  and  EUV  radia¬ 
tion.  The  level  and  intensity  of  such  radiation  is  directly 
proportional  to  the  level  of  solar  activity.  Observation  of 
sunspot  numbers  and  the  10.7  cm  solar  radio  flux,  shown  in 
Figure  10,  provide  a  means  of  predicting  ionospheric 
parameters,  such  as  the  F  region  critical  frequency  and  peak 
layer  height  [Smith  and  King,  1981] .  An  increase  in  solar 
activity  increases  the  horizontal  pressure  gradients,  which 
is  the  main  driver  of  the  thermospheric  winds.  At  the  same 
time  an  increase  in  the  charged  particle  population  leads  to 
increased  ion  drag,  which  inhibits  the  winds.  In  general, 
as  the  level  of  solar  activity  decreases,  daytime,  poleward 
wind  speeds  increase  because  ion  drag  decreases  more  rapidly 
than  the  pressure  gradient  force  [Rishbeth,  1972;  Babcock 
and  Evans,  1979;  Hedin  and  Mayr,  1987;  Lakshimi  et  al., 
1988]. 

In  Figure  11  and  Figure  12  the  dashed  lines  represent 
calculated  meridional  winds  for  a  2 4 -hour  period  based  on 
the  average  conditions  of  the  ten  most  quiet  days,  as 
recorded  in  Appendix  B.  The  solid  lines  represent  the 
statistical  error  as  calculated  at  3-hour  intervals.  The 


Fig.  10.  Variations  in  the  level  of  solar  activity  as 
seen  in  the  10.7  cm  solar  radio  flux  and  mean  sunspot 
number.  Units  of  flux  are  10'22  w  nf2  Hz'1. 
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decision  to  use  a  3-hour  interval  representation  is  based  on 
the  fact  that  this  investigation  is  primarily  geared  at 
studying  large-scale  cyclical  trends  of  the  meridional  winds 
as  opposed  to  looking  at  smaller  scale  phenomenon. 

The  results  of  this  investigation  indicate  that  the 
above  statements  are  essentially  correct.  Figure  11 
presents  meridional  winds  for  a  chain  of  ionosonde  stations 
which  stretches  from  43°S  to  62°N  geographic  latitude  and 
centered  roughly  at  140°E  geographic  longitude.  In  the  case 
of  the  Northern  Hemisphere  stations  of  Akita  and  Khabarovsk, 
the  poleward  wind  speeds  (positive  values)  increase  by 
roughly  30  m/s  from  solar  maximum  (1979)  to  solar  minimum 
(1987)  conditions.  It  is  noted,  however,  that  maximum 
poleward  wind  speeds  are  reached  roughly  three  to  four  years 
prior  to  solar  minimum  and  then  level  off.  In  the  case  of 
the  higher  latitude  station  of  Yakutsk,  an  increase  in 
poleward  wind  speeds  of  roughly  15  m/s  is  observed  between 
1979  and  1983.  A  gradual  decline  in  the  poleward  wind  then 
occurs.  This  result  appears  to  be  limited  to  this  station 
alone  and  may  be  due  in  part  to  variations  in  the 
calculation  of  hmax  by  the  IRI  model  ionosphere  and  that  of 
the  Bradley-Dudeney  model.  This  point  will  be  discussed  in 
more  detail  in  the  error  analysis  section  of  this  chapter. 

Figure  12  presents  meridional  winds  for  a  longitudinal 
chain  of  ionosonde  stations  centered  at  approximately  50°N 
geographic  latitude.  These  stations  also  exhibit  a  30  to  50 
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Fig.  11.  Variations  of  the  Meridional  wind^  for  a 
latitudinal  chain  of  ionosonde  stations  at  140T! 
longitude.  Cubed  lines  represent  calculated  vind 
velocities  and  solid  lines  represent  statistical 
uncertainties . 
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m/s  increase  in  the  speed  of  the  poleward  wind  as  time  moves 
from  solar  maximum  to  solar  minimum  conditions.  As  a 
general  rule,  these  stations  also  reach  a  maximum  poleward 
wind  speed  three  to  four  years  prior  to  the  time  of  solar 
minimum.  This  feature  corresponds  to  a  period  of  a  slower 
rate  of  decline  in  the  level  of  solar  activity  from  solar 
maximum  conditions  to  solar  minimum  conditions,  as  shown  in 
Figure  10. 

In  the  case  of  the  Southern  Hemisphere  stations  of 
Figure  11,  the  change  in  poleward  wind  speeds  is  not  as 
apparent  as  has  been  noted  for  the  northern  stations.  This 
is  particularly  true  of  the  station  of  Townsville  which  is 
closer  to  the  magnetic  equator  than  the  other  stations  of 
this  investigation  and  may  be  subject  to  additional  low 
latitude  effects.  Further  discussion  of  this  station  will 
be  deferred  to  the  section  on  the  repetitive  nature  of  the 
meridional  winds  relative  to  the  solar  cycle,  covering  the 
period  of  1965  to  1987.  Poleward  winds  for  the  most 
southerly  stations  of  Canberra  and  Hobart  again  demonstrate 
a  general  increase  in  wind  speed  with  decreasing  solar 
activity.  The  relative  increase  of  these  two  stations  is 
only  20  to  25  m/s  and,  as  previously  noted,  the  maximum 
speeds  are  reached  roughly  three  to  four  years  prior  to 
solar  minimum  at  1987. 

Corresponding  to  the  period  of  rapid  increase  in  solar 
activity,  1977  to  1979,  there  is  a  general  decline  in  pole- 
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ward  wind  speeds  of  roughly  30  to  35  m/s  for  the  Northern 
Hemisphere  stations.  A  more  drastic  decline  of  approximate¬ 
ly  40  m/s  is  noted  in  the  most  southern  stations  of  Canberra 
and  Hobart.  The  station  of  Townsville,  located  at  low 
latitudes,  is  again  the  exception  to  this  pattern. 

As  illustrated  in  Figure  11  and  Figure  12,  the  change 
in  the  equatorial  wind  speed  (negative  values)  of  the 
meridional  winds  show  a  more  drastic  solar  cycle  dependence 
than  the  poleward  winds.  An  overall  pattern  of  slower  wind 
speeds  occurring  at  solar  maximum  and  faster  wind  speeds  at 
solar  minimum  is  noted.  However,  the  slowest  wind  speeds  at 
all  stations  occurs  midway  through  the  solar  cycle  (1981  to 
1983).  In  the  case  of  the  Northern  Hemisphere  stations,  the 
winds  speeds  increase  over  a  range  of  70  to  160  m/s  between 
the  period  of  1983  to  that  of  solar  minimum,  at  1987.  For 
the  period  of  solar  maximum,  1979,  to  the  midway  point, 

1983,  the  meridional  winds  show  a  decrease  over  the  range  of 
40  to  130  m/s.  The  period  of  slowest  equatorial  wind  speeds 
corresponds  to  the  time  frame  where  the  slope  of  the  decline 
in  both  the  10.7  cm  solar  radio  flux  and  sunspot  number 
decrease,  as  seen  in  Figure  10.  The  Southern  Hemisphere 
stations,  with  the  exception  of  Canberra,  show  a  similar 
pattern  with  the  slowest  equatorial  winds  occurring  midway 
through  the  solar  cycle.  In  the  case  of  Hobart  and 
Townsville  a  range  of  increase  in  the  wind  speeds  of  70  to 
120  m/s  is  detected  for  the  period  of  1983  to  1987.  The 
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period  of  solar  maximum  to  midway  through  the  cycle, 
indicate  an  overall  decrease  over  the  range  of  100  to  150 
m/s.  The  station  of  Canberra  exhibits  a  difference  of 
approximately  25  m/s  between  both  periods  of  solar  maximum 
and  solar  minimum  versus  wind  speeds  calculated  at  the 
midway  period  of  the  solar  cycle. 

A  notable  difference  in  the  pattern  of  the  meridional 
winds  at  solar  maximum  versus  those  found  at  solar  minimum 
is  the  duration  of  the  equatorial  winds.  As  seen  in  Figure 
11,  this  variation  is  most  noticeable  in  the  Southern 
Hemisphere  stations  where  the  winds  blow  in  an  equatorial 
direction  for  approximately  15  to  20  hours  at  solar  maximum 
and  only  9  to  11  hours  at  solar  minimum.  This  result 
corresponds  to  the  work  of  Babcock  and  Evans  [1979]  for 
Millstone  Hill.  This  result  may  in  part  be  attributed  to 
the  difference  in  auroral  forcing  and  is  manifested  in  the 
summer  to  winter  interhemisphere  flow  regime  [Roble,  et  al., 
1977] . 

Repetitive  Efrtmre  af 

the  Meridional  Winds 

With  the  exception  of  Townsville,  all  other  stations 
exhibit  a  diurnal  fluctuation  in  the  critical  ionosonde 
frequency,  foF2,  and  in  the  associated  F2  peak  layer  height 
[Amaycenc,  1974;  Forbes  and  Garrett,  1976].  This  variation 
is  illustrated  in  Figure  13,  which  provides  a  representation 
of  the  variability  of  the  critical  ionosonde  frequency  foF2 


Fig.  13.  Representative  diurnal  variation  in  the 
critical  frequencies,  foF2  and  foE,  and  M(3000)F2 
throughout  the  solar  cycle.  Solid  lines  represent  foF2 
and  foE  and  dashed  lines  represent  M(3000)F2. 
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throughout  a  solar  cycle.  Comparison  of  the  critical  fre¬ 
quencies,  foF2  and  foE,  along  with  the  transmission  factor 
M(3000)F2,  and  the  solar  indices  for  the  cycle  of  1965  to 
1976  to  that  of  the  cycle  of  1977  to  1987  indicate  a  high 
degree  of  similarity.  Within  the  limits  of  uncertainty  of 
the  present  method  of  determining  the  magnitude  of  the 
meridional  winds,  one  would  expect  similar  patterns  and 
solar  cycle  variations  for  the  earlier  cycle  as  have  been 
detected  in  the  1977  to  1987  cycle. 

As  seen  in  Figure  14,  Townsville  exhibits  a  semi¬ 
diurnal  variation  in  the  ionosonde  parameter  of  foF2.  This 
variation  is  subsequently  reflected  in  the  values  of  hmax 
and  ultimately  in  the  magnitude  of  the  meridional  winds 
[Harper,  1973;  Hong  and  Lindzen,  1976;  Harper,  1979;  Mayr  et 
al.,  1979;  Burnside  et  al.,  1983].  Figure  15  illustrates 
the  semidiurnal  nature  of  the  meridional  winds  as  well  as 
their  repetitive  solar  cycle  variations  over  the  period  of 
two  solar  cycles.  As  has  been  discussed  previously,  both 
the  equatorward  and  poleward  wind  speeds  increase  with 
decreasing  solar  activity.  The  magnitude  of  these  varia¬ 
tions  is  less  than  for  those  stations  which  are  more  distant 
from  the  magnetic  equator.  This  is  most  likely  due  to  a 
need  to  account  for  the  E  x  B  zonal  wind  field,  wave  propa¬ 
gation  effects  from  below,  and  perhaps  due  to  the  existence 
of  a  more  complex  pressure  gradient  effect  [Harper  1973; 
Harper,  1979;  Mayr  et  al.,  1979;  Forbes,  1982].  Additional- 


Fig.  14.  Seal -diurnal  variation  of  the  critical  fre¬ 
quencies,  foF2  and  foE,  and  M(3000)F2  for  Townsville 
throughout  the  solar  cycle.  Solid  lines  represent  foF2 
and  foE  and  dashed  lines  represent  M(3000)F2. 
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ly,  the  slowest  equatorward  winds  coincide  with  the  periods 
of  a  change  in  the  slope  of  the  solar  flux  (1969  to  1971, 
and  1981  to  1983)  and  sunspot  numbers,  as  was  previouly 
discussed.  The  duration  of  the  equatorward  winds  solar 
maximum  also  exceeds  that  found  at  solar  minimum. 

Error  Analysis 

Statistical  uncertainties  as  described  in  Chapter  III 
are  accounted  for  in  the  prior  discussion  on  solar  cycle 
effects.  The  nighttime  equatorward  winds  exhibit  the  high¬ 
est  level  of  uncertainty  due  to  the  effect  of  the  a'2 
dependency,  as  seen  in  the  third  term  of  equation  (9) ,  com¬ 
bined  with  the  larger  difference  in  ho  and  hmax  for  this 
time  period. 

Global  wind  patterns  derived  from  ionosonde  critical 
frequencies  and  IRI  model  parameters  are  similar  in  nature. 
One  of  the  parameters  provided  by  the  IRI  model  is  hmF2,  as 
determined  by  an  empirical  fit  of  ionosonde  data  of  foE, 
foF2,  and  M(3000)F2.  These  are  the  same  parameters  used  in 
the  Bradley-Dudeney  model  employed  in  this  study  of  the 
meridional  winds.  Results  of  both  models  are  valid  only  at 
mid-latitude  regions.  The  results  of  this  study  agree  with 
the  findings  of  the  study  of  Miller  et  al.  [1989]  in  that 
differences  between  the  IRI  values  of  hmF2  and  the  Bradley- 
Dudeney  model  are  most  likely  the  result  of  local  phenomena 
in  terms  of  ionosonde  data.  Such  transient  behavior  is  too 
rapid  to  be  adequately  described  by  these  empirical  models. 
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Differences  in  these  two  models  is  most  noticeable  in  the 
calculation  of  nighttime  equatorial  winds.  These  differen¬ 
ces  do  not  change  the  findings  of  this  investigation  as  IRI 
calculated  nighttime  winds  generally  exceed  the  Bradley- 
Dudeney  values  at  solar  minimum.  This  would  tend  to  in¬ 
crease  the  magnitude  of  the  wind  speeds  at  this  time  of  the 
solar  cycle.  The  difference  in  hmF2  is  reversed  at  solar 
maximum  with  the  Bradley-Dudeney  values  of  hmax  exceeding 
the  IRI  model  values.  This  would  tend  to  reduce  the  calcu¬ 
lated  wind  speeds  and  lead  to  an  even  greater  variation  of 
the  meridional  winds  throughout  the  solar  cycle. 
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CHAPTER  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

At  mid-latitudes,  absorption  of  UV  and  EUV  radiation 
are  two  of  the  principle  ionization  and  heating  sources 
within  the  thermosphere.  The  level  of  such  radiation  is 
directly  related  to  the  amount  and  intensity  of  solar  activ¬ 
ity.  As  with  temperature  variations,  which  create  horizon¬ 
tal  pressure  gradients,  which  are  the  major  driving  force  of 
mid-latitude  thermospheric  winds,  there  is  a  definite  solar 
cycle  dependence  on  the  magnitude  and  direction  of  the 
meridional  neutral  winds.  Use  of  ionosonde  data  provides  a 
valuable  means  of  calculating  the  magnitude  and  direction  of 
the  global  scale  variations  of  the  these  winds. 

An  increase  in  solar  activity  increases  the  pressure 
gradients,  which  drive  thermospheric  winds,  but,  at  the  same 
time,  an  increase  in  the  charged  particle  population  leads 
to  an  increased  ion  drag,  which  inhibits  the  winds.  In 
general,  as  the  level  of  solar  activity  decreases,  daytime 
wind  speeds  increase  because  ion  drag  decreases  more  rapidly 
than  the  pressure  gradient  force.  One  area  of  future  study 
in  this  area  would  be  the  investigation  of  the  effect  of 
combined  auroral  forcing  and  solar  forcing  on  the  mid- 
latitude  thermospheric  winds.  This  type  of  study  would 
provide  a  much  needed  link  between  mid-  and  high  latitude 
ionospheric  models.  A  similar  study  could  also  be  conducted 
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at  low  latitudes  in  order  to  complete  a  truly  global 
ionospheric  model  of  thermospheric  circulation.  At  present 
the  computer  cost  and  time  requirements  may  be  prohibitive 
of  such  a  large  scale  study. 

For  the  slow  declining  phase  of  the  solar  cycle,  this 
study  indicates  that  for  the  Northern  Hemisphere,  as  the 
level  of  solar  activity  decreases  the  poleward  wind  speeds 
increase  by  30  to  50  m/s.  This  statement  is  illustrated  in 
Figure  16  which  provides  a  representation  of  the  peak  equa- 
torward  wind  and  poleward  winds  relative  to  the  10.7  cm 
solar  radio  flux.  As  a  general  rule,  the  poleward  winds 
reach  a  near  maximum  velocity  roughly  three  to  four  years 
prior  to  solar  minimum  and  then  level  off  until  solar  mini¬ 
mum  occurs.  This  result  varies,  to  some  degree,  from  sta¬ 
tion  to  station,  as  seen  in  the  case  of  Poitiers  as  illus¬ 
trated  in  Figure  16.  In  the  case  of  Southern  Hemisphere 
stations,  a  less  dramatic  increase  in  poleward  wind  speeds 
with  decreasing  solar  activity  is  observed.  For  the  South¬ 
ern  Hemisphere  locations,  results  of  this  study  indicate  a 
20  m/s  increase  in  the  poleward  wind  speed  is  achieved  three 
to  four  years  prior  to  solar  minimum. 

In  the  rapidly  increasing  phase  of  the  solar  cycle, 

1977  to  1979,  a  30  m/s  decrease  in  poleward  winds  occurs  in 
the  Northern  Hemisphere,  as  illustratec  in  Figure  16.  Mid- 
latitude  stations  in  the  Southern  Hemisphere  exhibit  a  40 
m/s  decrease  in  the  poleward  wind. 
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As  seen  in  Figure  16,  the  equatorial  winds  exhibit  the 
same  general  pattern,  as  seen  in  poleward  winds,  with  the 
faster  speeds  occurring  at  solar  minimum  and  slower  wind 
speeds  occurring  at  solar  maximum.  However,  a  variation  in 
this  pattern  is  noted  in  that  the  slowest  equatorial  wind 
speeds  generally  occur  midway  through  the  cycle.  This 
corresponds  to  a  decrease  in  the  rate  of  the  10.7  cm  solar 
radio  flux  and  sunspot  number.  Further  research  into  this 
particular  phenomenon  is  warranted.  Additionally,  in  the 
Southern  Hemisphere  the  duration  of  the  equatorial  winds 
also  exhibits  a  solar  cycle  dependence  in  that  the  duration 
of  the  equatorial  wind  increases  with  increasing  solar 
activity  due  to  increased  auroral  forcing.  Additional 
research  into  the  relative  importance  of  auroral  forcing 
versus  solar  forcing  of  the  mid-latitude  thermospheric  wind 
is  required  to  explain  this  circulation  phenomena  more 
completely . 

The  results  of  this  study  indicate  a  definite  trend  in 
the  winter  solstice,  meridional  winds  which  correlates  to 
the  11-year  solar  cycle.  In  the  case  of  poleward, 
meridional  winds,  a  direct  correlation  has  been  observed  as 
the  faster  wind  speeds  occur  at  solar  minimum  and  slower 
wind  speeds  occur  at  solar  maximum.  For  the  equatorial 
winds,  an  offset  is  observed  in  that  the  faster  wind  speeds 
occur  at  solar  maximum  and  minimum  while  the  slowest  wind 
speeds  occur  midway  through  the  cycle.  Such  results  have  a 
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Fig.  16.  Variation  of  the  equatorial  (nighttime) 
and  poleward  (daytime)  meridional  winds  for  one 
solar  cycle  relative  to  the  10.7  cm  solar  radio 
flux.  Units  of  flux  are  10"“  W  nf*  Hz"1. 
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bearing  on  mid-latitude  plasma  transports  phenomenon. 

Statistical  uncertainties  are  largest  in  the  nighttime t 
eguatorward  winds  due  to  the  effect  of  the  a'2  dependency  as 
seen  in  the  third  term  of  equation  (9)  of  Chapter  II.  The 
size  of  the  uncertainty  is  enhanced  by  the  greater 
difference  in  ho  and  hmax  for  this  time  period. 

Additional  areas  of  study  in  this  area  of  thermospheric 
circulation  might  include: 

1)  A  study  of  22-year  cycle  effects, 

2)  A  comparison  of  the  variability  of  thermospheric 
winds  with  solar  flare  activity  throughout  the  solar  cycle, 
and 

3)  Development  of  and  incorporation  of  a  wind  model 
for  the  current  mid-latitude  ionospheric  models  to  incorp¬ 
orate  ionosonde,  radar,  or  satellite  data,  which  are 
required  for  such  calculations. 
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Appendix  A.  General  Information 
on  Ionosonde  Stations 

Table  2.  Description  of  Ionosonde  Stations. 


NAME 

CEOCRAPHIC 

MAGNETIC 

L- SHELL 

LMT 

LOCATION 

LONGITUDE 

(hrs) 

(degrees) 

(degrees) 

Akita 

39.70  M/140.10  E 

206 

1.351 

BT+9.10 

Canberra 

35.30  S/149.00  E 

225 

1.955 

UT+10.  34 

Hobart 

42.90  S/147.30  E 

224 

2.629 

UT+10.34 

Khabarovak 

48.50  N/135.10  E 

200 

1.645 

UT+8.74 

Ottawa 

45.40  N/284.10  E 

351 

3.347 

UT-5.19 

Foitlera 

46.57  N/  0.35  E 

81 

2.408 

UT+0.83 

Sverdlovsk 

56.70  N/  61.10  E 

141 

2.357 

UT+4.77 

Tomsk 

56.50  N/  84.90  E 

159 

2.125 

UT46.03 

Townsville 

19.63  S/146. 85  E 

219 

1.321 

UT+10.00 

Yakutsk 

62.00  N/129.60  E 

194 

2.597 

UT+8.32 

LST 

MERIDIAN 

SUNRISE 

LMT/UT 

SUNSET 

LMT/UT 

NO  IONOSONDE 
DATA 

(degrees) 

135  E 

(bra) 

7.05/21.95 

(hrs) 

16.95/7.85 

(yra) 

None 

150  E 

5.14/18.80 

18.77/8.43 

1984  -  87 

150  E 

4.83/18.49 

19.20/8.86 

1984  -  87 

135  E 

7.39/22.65 

16.61/7.87 

1968,  76,  77. 
80,  82  -  83 

75  W 

7.20/12.39 

16.80/21.99 

None 

15  E 

7.20/  6.37 

16.80/15.97 

1966 

60  E 

7.88/  3.11 

16.00/11.23 

1982  -  83 

90  E 

7.88/  1.85 

16.17/10.14 

1972,  76,  79, 
81,  83 

150  E 

5.64/19.64 

18.42/8.42 

1983  -  87 

135  E 

8.36/  0.04 

15.64/7.32 

1965,  80  -  87 

Appendix  B.  Geomagnetic 
and  Solar  Data 


Table  3.  Geomagnetic  and  Solar  Data. 


YEAR 

F10.7 

F10.7A 

Ap 

Rl2 

1987 

69.6 

75.4 

4.5 

10 

1986 

72.9 

73.5 

4.3 

2 

1985 

72.6 

73.1 

6.3 

16 

1984 

100.3 

92.4 

5.2 

58 

1983 

140.8 

170.0 

7.3 

86 

1982 

152.2 

200.5 

4.2 

111 

1981 

161.6 

205.1 

4.8 

114 

1980 

201.4 

212.3 

4.4 

162 

1979 

189.2 

163.7 

7.8 

166 

1978 

94.7 

106.4 

3.5 

49 

1977 

74.1 

73.5 

5.2 

16 

1976 

70.9 

75.0 

5.6 

8 

1975 

75.3 

83.0 

5.4 

19 

1974 

86.8 

77.9 

5.5 

30 

1973 

97.6 

106.2 

5.6 

43 

1972 

101.1 

112.9 

4.7 

65 

1971 

150.0 

152.6 

3.5 

78 

1970 

148.9 

152.8 

4.2 

105 

1969 

158.9 

148.3 

2.6 

110 

1968 

189.2 

163.7 

6.7 

115 

1967 

161.6 

205.1 

3.0 

108 

1966 

94.7 

106.4 

2.0 

27 

1965 

74.1 

73.5 

2.6 

18 

16  February  1969 
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UNITED  STATES  DEPARTMENT  OF  COMMERCE 
National  Ocuwc  and  Atmospheric  Administration 
Environmental  Research  Laboratories 
325  Broadway 

Boulder.  Colorado  00303-3328 

E/CC2  D20354 
September  14,  1988:  DCS 


Dr  Kent  Miller 
Utah  Scare  University 
Center  for  Atmospheric 
and  Space  Sciences 
Logan,  Utah  84322-44-5 

Dear  Dr.  Miller: 

I  am  enclosing  the  final  part  of  your  data  re  st  (no.  E/CC2  D203S4) 

The  copies,  of  foF2,  FM3000F2 ,  and  foe  tab  sheets  are  all  that  are 
available.  Sow  stations  are  not  complete  and  you  nay  assume  the 
station  had  equipment  problems  and  was  not  operating  at  that  time. 

Please  be  advised  that  you  and  your  colleagues  are  always  welcome  to 
use  the  Boulder  Data  Center  when  you  are  in  town. 

Thnnlr  you  for  your  patience  In  getting  the  data  together  and  copied.  If 
you  find  any  of  It  not  readable  please  call  and  1  will  check  the  original 
data  to  see  If  we  can  help. 


Sincerely  yours, 

j.i) - 

Mrs.  Doris  Stansell 
Ionosphere  Branch 
NOAA/HCDC  E/GC2 
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